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THE FREQUENCY SHIFT OF A PULSE BY A TIME-INDEPENDENT, 
DISPERSIVE, LOSSY MEDIUM 


R.M. Jones 


ABSTRACT 


It is shown that a pulse undergoes a frequency shift while propa- 
gating through a medium whose losses depend on frequency. The amount 
of the frequency shift is proportional to the imaginary part of the 
complex group refractive index and the distance traveled in the mediun. 
No new frequencies are generated that were not in the original signal; 
the signal energy is merely "redistributed" by selective absorption. 
In an inhomogeneous medium, the frequency shift is proportional to the 
imaginary part of the complex group path. Furthermore, a beam propa- 
gating in a medium where the losses depend on the direction of propaga- 
tion will be bent toward the direction of least attenuation. Again, no 
new components in the angular spectrum are generated that were not in 


the original beam. The signal energy is merely redistributed by selec- 


tive absorption. 
1. SUMMARY 


When a wave pulse propagates through a medium in which the attenuation 
depends on the frequency, the pulse spectrum is attenuated faster on one side than 
on the other. The farther the pulse propagates into the medium, the more the 
frequency is shifted. 


This selective absorption might be interpreted as an effective Doppler 
shift of the pulse, even though the medium is time-independent. The effect 
depends partly on the time-dependence of the pulse; the effect vanishes for 


a monochromatic wave. 


The imaginary part of the complex group refractive index 


ar (nw) = i ((u-ix)w) 


is proportional to the freguency derivative of the absorption coefficient 


IY and is thus a quantitative measure of this effect. 
In an inhomogeneous medium, the imaginary part of the complex group path 
d 
' ez oe—— & 
P ry (w f n * ds) 


measures the cumulative effect of the frequency shift just as the real part 


of P' gives the time of travel of the pulse. 


Analogously, a beam that propagates in a medium that has angle-dependent 
losses will have its wave-normal direction shifted. This would be interpreted 
as an effective bending of the beam, even though the medium is homogeneous. The 
effect depends partly on the spatial inhomogeneity of the beam; the effect vanishes 


for a plane wave. 


As a particular example, the frequency shift of a time-harmonic wave modulated 
by a Gaussian-shaped envelope is inversely proportional to the square of the pulse 
length in addition to being proportional to the imaginary part of the complex group 


refractive index and the distance traveled through the medium. 


It is important to notice that all of the above effects are linear. That is, 
no new frequencies are generated that were not in the original signal. The signal 


energy is merely "redistributed" in frequency by selective absorption. 


John Bennett (1974, page 1583) announced these results as a private communi- 


cation from M. Jones. Here I document the results in full. 


2. BACKGROUND 


The propagation of pulses in lossless media is well understood (e.g., Stratton, 
1941; Hines, 195la; Panofsky and Phillips, 1955; Jeffreys and Jeffreys, 1956; 
Brillouin, 1960; Whitham, 1960; Budden, 1961; Jackson, 1962; Lighthill, 1965), in- 
cluding both isotropic and anisotropic media, in cases of anomalous dispersion 
where the group velocity no longer gives the velocity of energy flow, and the 
behavior of the precursor. For cases of normal dispersion, the velocity of energy 


flow is the same as the group velocity, 


v = 


" . (1) 


gg 


-> 
(Here, w is the radian frequency, r 21/A is the wave number, and k is normal 


to the phase fronts.) 


In general, a pulse changes shape as it propagates. There is, then, no single 
velocity that describes its motion. Under conditions of normal dispersion, the 
group velocity (1) gives the velocity of propagation of the pulse peak. However, 
the complete calculation of the propagation of a pulse reguires a calculation of 
the change in shape as well. The method is straightforward, although often diffi- 
cult in practice. A Fourier transform is performed on the initial pulse, the 
propagation of each Fourier component is found, and an inverse Fourier transform 
is taken to find the amplitude and phase of the signal as a function of time and 
position. Under conditions of normal dispersion, the inverse transform can often 
be evaluated by asymptotic methods (Jeffreys and Jeffreys, 1956; Price, 1968; 


Felsen, 1969). With care, the results are unambiguous. 


Many investigators have calculated the propagation of pulses of various 
shapes under various situations and calculated various aspects of the resulting 
distortion to the pulse (e.g., Felsen, 1969; Nicolis, 1967 (for a bell-shaped 
pulse); Wait, 1969a (exponential and bell-shaped pulse); Wait, 1969b (an example 
with a bell-shaped pulse that can be evaluated exactly); Vogler, 1969 (Gaussian- 
shaped pulse); Wait, 1970a (Gaussian-shaped pulse); Wait, 1970b (propagation of a 
Gaussian-shaped pulse over flat, imperfectly conducting earth); Wait, 1970c (Gaussian- 
shaped pulse); Vogler, 1970a, 1970b (Gaussian-shaped pulse)). 


The propagation of pulses in lossy media is less well understood. The main 
difficulty arises from trying to interpret the significance of the group velocity 
in (1) when it takes on complex values. Booker (1939) suggested that, since a 
pulse is made up of a spectrum of frequencies, a complex frequency might be found 
in the spectrum that would make the group velocity real. Hines (1951 a,b) showed 
that the velocity of the spatial maximum and the time maximum of a pulse do not 
usually coincide in a lossy medium. He argued that the velocity of the time maximum 


is more closely related to what is actually measured. He showed that 


S | (2) 


~ Teal(n') 


gives the velocity of the time maximum of a pulse, where 
(nw) = &— ((u-1x)0) (3) 
du) du) X 


is the complex group refractive index (Budden, 1961). (Here n is the complex 
phase refractive index [with real part u, and imaginary part - x], and c is the 
free-space speed of light.) 


Furutsu (1952), apparently unaware of Hines' work, argued that the concept 


of group velocity and wave path do not exist in a lossy medium. 


Suchy (1972 a,b,c, 1974) uses different arguments from Hines' to advocate 
real (30/ 8k) as a more appropriate group velocity. Suchy (1972a) argues that 
imaginary (30/ 9k) has no apparent physical meaning. Suchy (1974) interprets it 


in terms of a directional derivative of the real part of the wave number. 


The difficulty with interpreting a complex group velocity is similar to that 
of interpreting the complex ray direction that occurs for ray tracing in complex 
space (Poeverlein, 1962; Budden and Jull, 1964; Jones, 1970; Budden and Terry, 
1971; Keller and Streifer, 1971; Bertoni et al., 1971; Deschamps, 1972; Kratsov 
et al., 1974; Wang and Deschamps, 1974; Bennett, 1974; Connor and Felsen, 1974). 


The seemingly peculiar behavior of wave propagation in dissipative media has 


led various investigations into various aspects. For example, Hines (1951 c,d) 


and Arsaev and Kinber (1968) consider the direction of energy flux in dissipative 
media.  Poeverlein (1962) pointed out that absorption of waves can be represented 
by complex propagation vectors. Storey and Roehner (1970) and Roehner (1971) 
consider the direction of stationary phase for a beam of waves in an absorbing 
medium.  Bertoni et al., (1971) consider the nonlocal nature of propagation in 
lossy media.  Batorsky and Felsen (1971) consider complex waveguide modes.  Denman 


and Buch (1973) derive a Hamiltonian for dissipative systems. 


The practical calculation of the propagation of pulses in lossy media is 
not hindered by difficulties in interpreting a group velocity that takes on com- 
plex values. The method is the same as for lossless media, and the results for 
the amplitude and phase of the resulting signal as a function of time and position 


are just as unambiguous as for the lossless case. 


Several investigators have calculated various aspects of pulse distortion 
in lossy media.  Vogler (1969) calculated envelope shape distortion of a pulse 
propagating in a lossy troposphere.  Vogler (1970a and b) considered pulse dis- 
tortion further. Wait (1970b) considered pulse distortion for propagation over a 


flat, imperfectly-conducting earth. 


One particular aspect of the distortion of a pulse as it propagates through 
a lossy medium, that of the frequency shift of the carrier and the relation to 
the frequency dependence of losses and to the imaginary part of the complex group 
refractive index, does not seem to have received attention. John Bennett (1974), 
however, reported the result in the special issue of Proceedings IEEE on rays and 
beams, referring to some of my unpublished work. In the same issue, Connor and 
Felsen (1974), apparently independently, derive the same result for the case of a 
Gaussian-shaped pulse, but do not mention the relation to the imaginary part of 


the group refractive index. 


As we shall see, the above frequency shift is proportional to the product of 
the imaginary part of the group refractive index and the distance traveled in 
the medium. This discovery complements the already-known relation between the 
time of travel of a pulse and the product of the real part of the group refractive 


index and the distance traveled in the medium. 


Section 3 illustrates the main results in a gualitative way. Section 4 
establishes a quantitative relation between the imaginary part of n' and the fre- 
quency dependence of losses, and Section 5 derives a relation between the fre- 
quency dependence of losses and a frequency shift in the spectrum of a pulse. 
Section 6 applies the results to the case of a Gaussian-shaped pulse, Section 7 


generalizes to inhomogeneous media, and Section 8 discusses the propagation of 


beams. 


3. QUALITATIVE SUMMARY OF RESULTS 


Suppose a wave pulse 


Amplitude 


Time 


has a certain frequency and angular spectrum. 


PA 


Frequency or Angle 


Spectral Density 


When such a pulse propagates through a medium in which the attenuation is a function 


s 


Frequency or Angle 


of frequency and/or angle 


Attenuation 


the pulse spectrum will get attenuated (or "worn away") faster on one side than 


the other. 


Amplitude 


ti 
t2 
t3 
WA 


Frequency or Angle 


The result of frequency-dependent losses is that the farther the pulse pro- 


pagates into the medium, the more the frequency gets shifted. 


Amplitude 


Time 


This selective absorption might be interpreted as an effective Doppler shift 
of the pulse, even though the medium is time-independent. The effect depends 


partly on the time dependence of the pulse; the effect vanishes for a monochromatic 


wave. 


Such a "single pulse" Doppler shift has nothing to do with the Doppler shift 
measured by the pulse-to-pulse phase shift in a Doppler radar. This effect might 


influence range measurements in an FM-CW radar, however. 


The result of angle-dependent losses is that the farther the wave packet 


propagates into the medium, the more the wave-normal direction gets shifted. 


This could be interpreted as an effective bending of the wave packet, even 
though the medium is homogeneous. The effect depends partly on the spatial in- 


homogeneity of the wave packet; the effect vanishes for a plane wave. 


4. THE RELATION BETWEEN THE IMAGINARY PART OF n' AND THE 
FREQUENCY DEPENDENCE OF LOSSES 


Suppose a monochromatic plane wave with frequency w travels in the z direction 
through a homogeneous medium that has a complex phase refractive index n. The 


time and space variation of the amplitude A of the wave is 
id A : . Un 
Ae' = A, exp(id - az) = A exp(iut - A z) (4) 
where $ is the phase of the wave and a is the absorption coefficient. 


To obtain an explicit expression for the freguency dependence of the phase 


and amplitude, let us take the derivative of the phase with respect to freguency, 


do . = d = > | 
Fam t real x (nw) po real (n'), (5) 


and the logarithmic derivative of the amplitude with respect to frequency, 


a AP 43 2. d. MS Nr ' 
> (Log A) z = — imag 7, (nw) = 7 imag (8-2, (6) 


where we have used (3) for n'. Equation (5) leads to the well-known result (2) for 
the group velocity, while (6) gives the new result, 


_ da „ima n' 


dw c T 


(7) 


This means that the imaginary part of the complex group refractive index is pro- 
portional to the frequency derivative of the absorption coefficient. Evidently, 
the group refractive index is real if and only if the absorption coefficient does 


not depend on frequency (at least at a given frequency). 


5. THE FREQUENCY SHIFT. OF A PULSE 


Suppose we have a time-harmonic wave of frequency wo modulated by an envelope 


having one peak at t = 0: 


m(t) exp (iw t) = k f exp(iwt) A(w) exp(B(w))dw (8) 
2T —© 
where 
A(w) exp(B(w)) = — f  exp(-iwt) m(t) exp(iw_t)dt (9) 
: V2T =% ” 


is the Fourier transform of m(t) exp (iu t) and has been divided into a slowly 
varying part A(w) and a quickly varying part exp(B(w)) to show its asymptotic 
behavior. Let the pulse propagate through a homogeneous medium that multiplies 


a monochromatic wave by 
exp(- i 2 uz) exp (-az) (10) 
where u is the real part of the phase refractive index, G is the absorption coeffi- 


cient, and z is the distance traveled (along the wave normal direction) through the 


medium. The spectrum of the pulse at a point z in the medium is then 


A(w) exp(B(u) - i  u(u) z - aw) z) . (11) 


To find the dominant frequency u in thís spectrum, we set 


$- law) exp(B() - alw) 2|, -0- (12) 
| P 


When the slowly varying term A(w) is neglected, (12) gives 


d da (Ww ) 
real au, Py == a * (13) 


We now expand B(w) about w = W, > where w, is the dominant frequency of the trans- 


1 
mitted pulse. That is, we expand about the point where 


d B(w) 
real ———— = real B'(w.)=0. | (14) 
du, 4 


(Notice that if B is an analytic function, the imaginary of the above is also zero, 


so that w is then a saddlepoint of the transmitted pulse.) (In many cases, 04 = W+) 
The Taylor expansion gives 
B(w) = B(w ) + B' (w. )(w-w,) + Žau ) (0-0 j^ + (15) 
1 1 1 2 1 i pater: 


where primes indicate derivatives with respect to the argument. Substituting (15) 
into (13) gives 
d a(w_) 
" D M E | 
real B (0, ) (v, 04) Ee ura 0. (16) 
p 
We have neglected second and higher powers of (w - w) to give the asymptotic 
high-frequency behavior. Solving (16) for the frequency shift, a ~ 04); gives 


: d a(w ) imag n'(w ) 
P + real B (w) d A c real B (w) 


We have used (7) to combine with the result of the previous section. Thus, the 


medium shifts the frequency of the pulse by an amount proportional to the distance 


10 


traveled in the medium and proportional to the freguency derivative of the ab- 
sorption coefficient (and thus proportional to the imaginary part of the complex 
group refractive index). (It is probably accurate enough in most cases to use © 


n' (w,) for n' (w) ia (17).) 


The physical explanation for this frequency shift is clear. If the absorption 
depends on frequency then the medium will selectively absorb either the higher or 
the lower frequencies in the pulse spectrum, thus shifting the dominant frequency 
in the spectrum. It is also clear that thís frequency shift should be proportional 
to z, the distance traveled in the medium. It is also reasonable that the frequency 
shift is larger, the larger the spectral width of the pulse, which is inversely 


proportional to the real part of B", the curvature of the spectrum peak. 
6. EXAMPLE -- A GAUSSIAN-SHAPED PULSE 
For a Gaussian-shaped pulse, let 
r 
m(t) = exp(- (t/1)^) (18) 


where T is a measure of the pulse length. 


(Although the Gaussian-shaped pulse is not realistic because it has no beginning, 
the propagation of the peak of a Gaussian-shaped pulse should be similar to that 


of realistic pulses with rounded peaks.) 


Then from (9), the pulse shape in the frequency domain is 


2 
Aw) exp(B(u)) = — expC- j- (99 )^) , (19) 
2 
and thus 
B(w) = Es é 
U > 7 0 (20) 
and 
W, 7 0. . : (21) 


41 


Substituting the above into (17) gives the freguency shift of the spectral peak 


imag n'(u ) | 
W =W ee | (22) 


Equation (22) shows that the frequency shift is larger for shorter pulse lengths 


(i.e., a wider pulse spectrum), as we expected. 


[ INHOMOGENEOUS MEDIA 
For inhomogeneous media, (10) can be replaced by the transfer function 
a(w) exp(- i > P(w)) ; (23) 


where a is a slowly varying function of w, and P is the phase path or phase integral 
(Budden, 1961, page 136). Thus, neglecting the slowly varying term a(w), (16) 
would be replaced by 


real B" (03) (05704) += imag P! (au) -0, (24) 
where 
P'(u) = $- Ww PW)) (25) 


is the group path or equivalent path (Budden, 1961, page 280). Solving (24) gives 


the frequency shift 


. imag P'(u ) 


W.-W : | (26) 


DI Te real B" (w) 


Equation (26) is a general relation that applies to propagation through any medium 


that has a transfer function of the form (23). 


For example, it would apply in the ray-theory approximation to waves propagat- 


ing through an inhomogeneous anisotropic medium if we use 


24 


Pa Pond. | (27) 
ray 
path 
P can be called the generalized or complex phase path (Budden, 1961, pages 173, 
197, 507), n is a vector pointing in the wave-normal direction and having a value 
equal to the complex phase refractive index (Budden, 1961, page 531), s is the 
distance along the ray path, and ds is a vector pointing in the ray direction. 


As another example, a full-wave reflection or transmission coefficient could 
be put into the form of (23). Then applying (26) would give the frequency shift of 


a pulse having a spectrum given by (9). 


8. BENDING OF A BEAM IN A LOSSY HOMOGENEOUS MEDIUM | 


The results just presented show that a pulse undergoes a frequency shift 
while propagating through a medium whose losses depend on frequency. We ordinarily 
expect frequency shifts only in time-varying media, so obtaining a frequency shift 


in a stationary medium is surprising. Two conditions are necessary for this result: 
1) A wave has a finite frequency spectrum. 


2) The medium has freguency-dependent losses. 


It is interesting to note that no new freguencies are generated that were not 
in the original signal. The signal energy is merely "redistributed" by selective 


absorption. 


It is tempting to ask whether these results, derived for the freguency 
spectrum of a pulse, apply to the angular spectrum of a beam. More explicitly, 
will a beam bend in a homogeneous medium having losses that depend on the propaga- 
tion direction of the beam? The improbable answer is yes, because components 
of the beam on one side of the angular spectrum will be attenuated faster than 
those on the other side as the beam travels through the medium, giving a curvature 
to the beam. However, the beam does not actually bend to directions where there 


was no wave energy to begin with. 


13 


Thus, in analogy to the above results, if the group velocity 


12 


(28) 


ele 


is complex, then the medium has losses that depend on the direction of propagation, 


and a beam propagated in such a homogeneous medium will be bent toward the direction 


of lowest attenuation. 
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